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Abstract 

Two of the approaches to the hadronic productions of the double heavy 
mesons Be and B* are investigated. Comparison in various aspects on the 
results obtained by the approaches is made and shown in figures and a table. 
Some trial understanding of the approaches themselves and the achieved re- 
sults is presented. The results may be used as some references for discovering 
the mesons at Tevatron and LHC. 
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Recently the interest in the Be meson, one of the double heavy flavor mesons, is 
aroused widely due to its properties. Similar to the heavy quarkonia, rjc, J/ip, and 
rjb, T etc., it is a double heavy quark- antiquark bound state, so the QCD inspired 
potential model will work well for describing the binding effects of it'^'^l ; but different 
from them, it carries flavors explicitly, so it may decay by weak interaction only, and 
as a result, it has a much longer lifetime ( a typical weak decay one ) and plentiful 
decay channels which have sizable branching ratios Especially, some of its 

decays can be calculated quite reliably and may be measurable in the near future 
Thus the meson in addition to the heavy quarkonia, may be used to test the 
QCD inspired potential models and the acting weak decay mechanisms for relevant 
heavy flavors further. Another important reason to make the B^ physics interesting, 
is the study of the Be meson being accessible soon experimentally. As pointed out by 
several independent theoretical estimates'^"^^', the cross sections of its production 
at certain existent and planned colliders are sizable and some typical signals may 
project over the background. 

Having all the possible productions of the double heavy flavored meson reviewed, 
the authors of refs.[5,6] have pointed out that the most suitable ones of high energy 
processes to produce sufficient events of the Be meson at the existent and planned 
facihties, are those at a Z° boson 'factory', such as LEP-I, and of energetic hadronic 
collisions at Tevatron and LHC etc. In ref.[5], besides a complete calculation on 
the Be meson production at the level of the lowest order of perturbative QCD 
(pQCD), the fragmentation functions for b ^ Be and b B*(S-'wme) were also 
worked out correctly, while those for b — > X(bc) (-P-wave) in ref.[7]. In fact, it is 
the first time to work out the fragmentation functions correctly, because not all the 
terms, being the lowest order, had been taken into account until the authors did. 
The fragmentation functions obtained by ref.[5] were confirmed by others soon'^'^^1 
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i.e. the authors of refs.[8,ll] recalculated the same fragmentation functions out 
in an axial gauge, a different gauge from adopted in ref.[5], and the factorization 
of the fragmentation functions are manipulated more manifestly. Furthermore the 
evolution of the fragmentation functions with changing of the fragmentation energy 
scale was also considered in refs.[8,ll] by solving the corresponding Altarelli-Parisi 
equation'^^l . 

Of all the proposed theoretical estimates^^"-'^^! , besides those of pure phenomeno- 
logical ones with Monte Carlo simulation such as done in ref.[10], the adopted ap- 
proaches for estimating the hadronic productions of the Be meson may be divided 
into two categories, although they produce results in consistency in order of mag- 
nitude (different in values from each other), and all are based on pQCD. The first 
category is to consider the production in a fragmentation picture i.e. the Be meson 
is produced due to fragmentation of a heavy flavor jet (here b jet mainly) 1^1. It is 
very similar to that of a light meson production from a jet, but the fragmentation 
energy scale is much higher (above that of nonperturbative QCD) that the fragmen- 
tation functions are calculated with pQCD. According to pQCD, the production 
cross section: 

ijk ^-[^■^ 

■daij^kx{xiX2X3, hf) ■ D^'{x3,hf), 

FIj{x, hf) is the distribution function of the parton i in the hadron i7, ddij^kxi,- • ■) 
is the cross section for the relevant jet inclusive production {i + j ^ k -\- X) and 
D^^ix^jiF) is the fragmentation function of Be from jet k. The formulation here 
means that the calculation should be carried out at a typical energy scale [If of the 
process. The specific fragmentation functions D^''{x, /if), {k = b, c) are calculated in 
the framework of pQCD. In this approach, it is easy to extend straightforwardly up 
to the leading logarithm approximation (LLA) accuracy level. The second category 
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is not, as done in the first category, to factorize the "subprocess" i + j ^ Be + X 
into two factors further: the 'jet production' i + j ^ k + X and the 'fragmentation' 
of the Be meson from the jet k ^ Bc + X, but the subprocess is treated as a whole, 
and to compute it directly in the framework of pQCD too i.e. the production cross 
section: 

d(7 = Y^j^^lj dX2F}j^{xi,llF)Fjj^{x2,I^F)daij^BcX{XlX2,I^F)- (2) 
ijk 

Although the computations of the second category are available only up to the 
lowest order of pQCD so far, in principle, they may be extended to higher orders 
with lengthy and boring calculations. In both of the two approaches (at the lowest 
order approximation), the wave function at original of the (c6) bound state system 
will occur in the fragmentation functions and the amplitude of the subprocess i+j —>■ 
Bc+X respectively, whereas the wave function may be obtained from potential model 
for the double heavy bound state system precisely'^'^'^'^-'^l . As B*{l^Si) meson has 
a cross section for hadronic productions bigger than that of B^, and it will decay 
to the ground state Be with a branching ratio almost 100% in a very 'short' time 
(without decay vertex in experimental detector) so it contributes the B^ production 
substantially, thus in the paper we will discuss B* and Be together from now on. 

In hadronic productions of Be{B*), the substantial contribution is from the sub- 
process of gluon-gluon fusion g + g ^ Be{B*) + b + c but not from quark-antiquark 
annihilation q + q ^ Be{B*) + 6 + c at a relatively high energetic colliders such 
as Tevatron and LHC etc'^^, thus we will restrict ourselves to discuss the subpro- 
cess of gluon-gluon fusion and to find out the differences attributed to the adopted 
approaches from now on in the paper. For convenience, we will denote the first 
category as Approach-I, whereas the second one as Approach-II. 

Prom the knowledge of pQCD, Approach-I depends more on the factorization 
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theorem, whereas Approach-II, being much more comphcated than Approach-I, is of 
a fixed order complete calculation. If Approach-I is extended up to the level of LLA, 
it may achieve better results for very high energy and high Pt problems. Whereas for 
some of the other problems with not very high energy and very high Pt, the complete 
fixed order approach, Approach-II even at the lowest order may achieve better ones. 
We are interested in examining the two categories of the approaches quantitatively 
for the Bc{B*) production not only due to the the experimental interest, but also 
due to the theoretical interest, because the quantitative results may offer references 
for study of Bc{B*) experimentally and understanding the production mechanics 
theoretically. 

The precise differences between the two categories in the production of the dou- 
ble heavy meson Bc{B*) in decay, have been given in ref.[5], although there the 
LLA corrections for Approach-I were not considered and the results for the rate of 
Bc{B*) production were underestimated owing to a smaller QCD coupling {as{Q'^) 

= "^|o) being adopted. Since the relevant 6-jet produced in decay is very 
energetic {Ei, — at C.M.S.) and the process is comparatively simple, thus 

to find the correspondece of the approaches is simple, the difference in values be- 
tween the approaches in partial width is less than 20%. However, in the hadronic 
productions the colliding energy of the subprocess varies in a wide region and may 
be quite 'low' at Tevatron, even at LHC, and the subprocess itselfis much more 
comphcated. The energies of the subprocess in the hadronic collision at Tevantron, 
even at LHC, in most chances are much smaller than that in Z° decay, and we will 
return this point more precisely later on. Furthermore the subprocess is much more 
comphcated than that of Z° decay: there are 3 diagrams instead of one of Z° decay 
for Approach-I, and there are 36 diagrams instead of 2 for Approach-II. Concerning 
the hadronic Bc{B*) production, the masses of the heavy quarks inside the meson 



5 



Bc{B*) play the role as the proper energy scale in the production, and it is much 
greater than IVqcd so pQCD calculations are always applicable, no matter how big 
Pt the heavy flavor (here the Bc{B*) meson) carries, that is very different from hght 
flavor productions. Generally to know the accurate contributions from the low Pt 
components quantitatively, which is expected to have substantial different results 
for the approaches, is interesting, because it is necessary for writing an event gen- 
erator which may produce reliable low Pt events. To see how low Pt events could 
be well detected by the concrete detector is interesting for the experiments on the 
concerned subject(s) i.e. Bc{B*) mesons for present problem, whereas without the 
reliable event generator, the aim could not be reached at. 

The differences between the two categories of the approaches, in fact, may be 
attributed how to deal with the subprocess of Approach-II. In Approach-II we deal 
with it as a whole i.e. a complete pQCD calculation on the process, though only the 
lowest order one is available so far; whereas in Approach-I it is treated to produce 
heavy quark jets flrst and then to fragment a meson Bc{B*) from one of the pro- 
duced heavy jets, thus as known from the proof of the pQCD factorization theorem, 
Approach-I is not as good as Approach-II even doubtable, if the jet responsible for 
the fragmentation of a Bc{B*) meson, is not very energetic. 

To understand the differences of the approaches, let us analysize the subprocess, 
g + g —>■ Bc{B*) + b + c, carefully. According to Approach-II, there are 36 Feymann 
diagrams responsible for it. Some of the typical ones are collected in Fig.l(a,b). It 
is easy to realize that the diagrams (the amplitude) may be divided (splitted) into 5 
independent subgroups (terms) according to their color structure, and each of them 
alone is guage invariant '^l It is too long to write down here the total amphtude of 
the subprocess explicitly, however, we may write its color structure out explicitly in 
a short formulation, and with it we will be able to find out some correspondence 
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and difference of the two approaches. In general, the formulation for the amplitude: 

6 

A{a,b,t,j) = ^C:^,.M„(ei, £2,^1,52). (3) 



a=l 



Here each of the color factors C^^j{a — 1, 2, • • • , 6) is a product of the Gell-Mann 
matrices: 

— — (A" • X'')ij ; 

— — (A'' • A")jj ; 

—(A" • X'')ij ; 

— (A • A")jj ; 

<5,,tr(A«-A'')-l(A'^A%; 
1 
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However, we should note here that not all these color factors are independent, be- 
cause there exists a relation among them, that is 

/^ab /^ab /^ab /^ab /r\ 

Therefore only 5 color factors are independent, and we may choose them as: 

C'^% = C^ij (when m = l,---,4); 



r^'ab ryab r^ab 

*-^5 ij — *-^3 ij *-^5 ij- 



(6) 



Thus the amplitude may be rewritten as: 

Aia, b, i, j) = ^ C^?,M^(ei, £2, si, ^2). (7) 



Being independent, the coefficients of the color factor C'^lj, the sub-amplitudes 
Mj(. ( k — 1,2, •••,5 ) are individually gauge invariant, thus each of them may 



acquire certain meaning. Owing to the fact that each of the amphtudes is re- 
lated to certain Feynman diagrams of the 36 precisely, the explicit formulas of 
{ k = 1,2, ■ ■ ■ ,5 ) may be written down directly, based on the rules of the duel am- 
plitude method^^l Therefore one may find out the correspondences and difference 
between Approach-I and Approach-II. Since the Approach-I is of a fragmentation of 
a 6-quark jet in the Bc{B*) production (the fragmentation of a c-jet contributes too, 
but it is much less important than that of a &-jet), with the decomposition eqs.(4-7) 
and according to the color structure to trace back to the diagrams, one may find 
the correspondence: the sub-amplitudes C'^^^Ml (A; = 1, 2) is to correspond to some 
of Approach-I's amplitudes, whereas those of C'^^jM^ {k = 3,4,5) cannot find any 
correspondence in Approach-I. The fact of the correspondences is easy to be under- 
stood by means of the Feynmann diagrams of the two approaches: One may find 
that the diagrams such as Fig. 1(a) which contribute to the sub-amplitudes C^^f^M^ 
{k = 1,2) substantially in the sense of the color structure, could be understood 
as if two jets were produced and the fragmentation of Bc{B*) meson was followed, 
whereas for the diagrams such as Fig. 1(b), which contribute to the sub- amplitudes 
C'^^jMl {k = 3,4,5) substantially, there is no similar correspondence at all in the 
above sense to Approach-I. Therefore we expect the results achieved by the two 
approaches being different, so a thorough investigation of the approaches quanti- 
tatively, even though numerically, is interesting. We will devote this paper to the 
investigation]^, i.e. to compare the Bc{B*) hadronic productions of the approaches 
quantitatively in various aspects. We will plot the numerical results of each ob- 
servable, obtained by the two approaches into one figure together, different figures 

show different aspects of the approaches, and finally we will try to reach at some 
^ During the period of revising the paper, several papersl^^'^"' come out and certain disagree- 
ments on the calculations are presented, thus to clarify the situation is also necessary. 
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conclusions. 

First of all, we calculate the total cross sections of the Bc{B*) productions by the 
two approaches at various hadronic coUiders i.e. for various C.M.S's energies of the 
coUiding hadrons, but only the lowest order for the subprocess g+g Bc{B*) +6+c, 
is concerned. The obtained total cross sections are put into Tab.l. We should note 
here that throughout the paper without special statement, the following manipula- 
tions and parameters are taken. When calculating the productions of p — p and p — p 
collisions, only gluon-gluon fusion mechanism is considered due to its domination 
over the others such as quark-antiquark etc'^l; the CTEQ3 structure functions with 
A^^iuf = 4) = 0.239GeV (corresponding a,(m|) = 0, 112)[^41 and a^iQ^) with an 
energy scale = s/4 (s is the cm. energy squared of the subprocess) are adopted. 
As for the masses, the values nic — l.bGeV, rub — 4:.9GeV and Mb^(b*) — QAGeV^'^^ 
are taken. Furthermore in the calculations the wave functions of Be and B* at origin 
are obtained obtained from potential model and the difference of wave functions, as 
the masses, for Be and Bq is ignored here for the 'lowest order calculation'. In order 
to compare with those adopted in literature easy, what we adopt it here is in decay 
constant formulation: /^^ ~ 480MeV (under the convention /^r = 132MeV). 

Note here that in the table when calculating the subprocess g+g Bc{B*)+b+c 
at VI — 20, 30, 60GeV, a constant of strong coupling q;^ — 0.2 is taken, and when 
the row is denoted with a '*' ('**')' the results are indicated to have a cut for small 
Pt < 5GeV {Pt < lOGeV). Up to the concerned order of pQCD, the uncertainties 
here come only from the choices of the values of rric, nib, Mb^{b*), cks and /b^. 

TABLE I. The total cross sections for the productions of the Be meson 
and its excited state B* obtained by the two approaches (in unit nb). 
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Collision 


Approacli-I 


Approach-II 




(lifi-o) Bl {l^Si) 


B, (li5o) Bl (l35i) 


pp{VS= l.STeV) 


0.747(4) 1.23(1) 


0.850(8) 2.07(2) 


ppi^/S = l.STey,*) 


0.229(2) 0.389(3) 


0.259(4) 0.646(6) 


pp{y/S =1.8TeV,**) 


0.0331(9) 0.0570(6) 


0.0373(1) 0.0894(3) 


pp{s/S = lATeV) 


8.63(5) 14.0(1) 


10.6(1) 26.4(3) 


pp(\/5 = 14TeV,*) 


3.07(3) 5.11(4) 


3.71(6) 9.43(9) 


pp{sfS= MTeV,**) 


0.584(7) 0.986(10) 


0.698(1) 1.69(4) 


gg{^f^ = 2{iGeV) 


0.704(5) • 10-2 0.118(1) • 10-1 


0.661(7) • 10-2 0.160(2) • 10-^ 


gg{V^ = SOGeF) 


0.678(8) • 10-2 0.103(1) • 10-^ 


0.949(8) • 10-2 0.244(3) • 10-^ 


gg{^s = eOGeV) 


0.321(7)- 10-2 0.456(9) • 10-2 


0.782(9) • 10-2 0.203(3) • 10-^ 



The Pt dependence of the productions at various coUiders Tevatron and LHC 
is interesting experimentally, thus we have calculated it and plotted the results 
in Fig. 2. In the calculations, the low Pt component contribution has been taken 
into account too, though for Approach-I the computation is problematic. It is 
because the production closing to the threshold (where Pt cannot be big) needs 
special consideration and corrections in Approach-I, but here we merely make the 
'approximation': the Pt of Bc{B*) being fixed in the direction of the produced 
heavy quark jet, in fact, it is not a good approximation when the 'fragmentation' 
is very close to the threshold of the Bc{B*) meson production, thus the low Pt 
component contribution as shown in Fig. 2 is not so well estimated for Approach-I. 
From the figure one may see that for the Be production, the difference between the 
two approaches is not sizable but for the B* production it is quite great (about a 
factor two even greater) and, general speaking, as Pt is going high the production 
cross sections predicted by the two approaches are approaching to equal (for B*, up 
to Pt = 20GeV they are still different). 

In order to have an outline about the gluon-gluon subprocess in hadronic coUi- 
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sions, in Fig. 3, we present the production cross sections at Tevatron and LHC versus 
the colhsion energy s of the glouns inside the colhsion hadrons. As the small Pt 
component of the productions is not able to measure, we have imposed a cut for 
those of small Pt{^ 5GeV) here. From the figure, one may see the cross sections 
drop in a logarithm scale versus s increasing. In fact, if we had not imposed the cut 
for small Pt, the cross sections would have a "peak" around 20GeV ( not very far 
from the threshold of the subprocess ~ 12.8GeV ), and then would drop 0. One 
may see that when s reaches at 80GeV, the cross sections have dropped down at 
least one more orders of the magnitude already. Thus one may conclude that in the 
hadronic collisions the dominant contribution to the Be and B* meson productions 
is not from very high energetic gluon fusion but from relatively low energy, that is 
the great difference from that in decay as we emphasized earlier in the paper. For 
B* production, the cross sections obtained by Approach-II are greater than those 
by Approach-I at various energies with a factor 5 or greater, but for B^ production, 
the difference caused by the two approaches is within a factor 2, less than that for 
B* production. 

We should note here that besides the cut for small Pt being imposed and the cou- 
pling constant as being running, the cross sections in Fig. 3 are different in meaning 
from that of the gluon-gluon fusion for precise a/I in Tab.l, as the later is merely of 
gluon-gluon fusion but the former has the structure functions of the collision hadrons 
convoluted into. 

All the resultant cross sections of the hadronic productions are achieved always 
by a convolution of the cross section of the relevant subprocess and a common factor, 

the structure functions of the incoming hadrons of the collisions. In order to highlight 

■^To shorten the paper and to present the more useful results, we would not present the curves 
without Pt cut here, although we have them. 
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the differences of the two approaches, we have also calculated the subprocess cross 
sections as if the subprocess is an independent one, i.e. the cross sections of the 
gluon-gluon fusion at various precise energies. The total cross sections have been 
put in Tab.l already, but the transverse momentum Pt and rapidity Y distributions 
at various C.M.S's energies are presented in Figs. 4: in Fig.4(a,b) for = 30 GeV, 
Fig.4(c,d) for = 60 GeV respectively. In these calculations, we have taken 
as = 0.2, a constant, as emphasized above. One may see the fact very clearly that, 
as expected, the values obtained by Approach-II at small Pt and small Y are always 
greater than those obtained by Approach-I, whereas the values are approaching close 
when Pt or Y increases. 

In summary, the two approaches cause some substantial differences in total cross 
sections and the Pt distribution etc. indeed, especially at low Pt- Approach-II 
should be suitable at low Pt, even at low Pt and low Y both. The heavy masses of 
the quarks inside the meson play the role to offer the least and proper energy scale in 
the concerned productions and to guarantee pQCD being applicable always. Namely 
it is the heavy quark masses being the least proper energy scale, instead of Aqcd, 
that appear in the formulas (appear in as, the coupling constant for the lowest order 
calculations and in the logarithmically large terms if higher order calculations are 
carried out). Furthermore, the most important productions of the Bc{B*) meson are 
shown in Fig. 3 not due to very energetic parton collisions, and from Fig. 2 one may 
also see that the Pt cannot be great in the interesting processes concerned in the 
paper, thus in order to collect as many as possible events of the mesons and B* 
so as to discover them and to study their properties, one could not estimate the low 
Pt components of the productions too rughly from very beginning and should try 
to have a good one as one can. For this purpose, it is sure that Approach-II is good, 
and the logarithmical terms to the heavy quark masses need not to worry about 
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too much, as the terms appear from high order calculations and become important 
at very large Pt only^^^l. Approach-I is better than Approach-II for estimating the 
productions at very large Pt if the former taking into account the large logarithmical 
terms to the heavy quark masses by LLA but the later not (as the present case). 
Whereas for the concerning hadronic productions at Tevatron, even at LHC, the 
former's advantages have not "matured" yet because of the same reason as pointed 
out above: Pt cannot be very great at the concerned processes. 

We think that the defferences of the obtained results by the two approaches may 
be understood by the fact that Approach-II involves more mechanisms than those 
Approach-I does, as argued in terms of eqs. (3-7) early. Recently the authors of 
ref.[17] also recognized that certain higher order gluon fragmentation besides the 
fragmentation of a heavy quark may contribute to Be production substantially. 

According to the experiences of heavy quark productions in hadronic collisions 
and the theoretical loop calculations, we know that a full perturbative QCD calcu- 
lation up to one loop level may achieve quite high accuracy'^^' , thus a higher order 
full perturbative QCD calculation on the hadronic production of the double heavy 
flavor meson Bc{B*) under Approach-II will be very interesting surely^^^l. 

In the procedure of revising the paper, several calculations'^^"^^' on the same 
problems appear. The authors of ref.[20] have found that in their earlier version 
they had omitted a color factor l/\/3 in amplitude, and when having the mistake 
corrected they have found a nice agreement between theirs and those of ref.[6]. 
The authors of ref.[19] have investigated various calculations quite systematically, 
therefore we have checked the numerical results for the subprocess g + g ^ Bc{B*) + 
6 + c at VI = 20, 406*6^, by means of our program but having their parameters. 
As a result, we have found that our results and theirs are consistent with each 
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other exactly within the Monte Calor errors^, that in fact is a confirmation for our 
programs and those of ref.[19]. Only for updating, the new version of the structure 
functions CTEQS'^^^ and the parameters appearing in the calculations as quoted 
above (e.g. the energy of HLC ^/S = UGeV etc), which are slight different from 
those of ref. [19] even ref.[6], have been adopted here so that the numerical results for 
the total cross sections and the other obsevables involving two structure functions 
of the two colliding hadrons are reasonably different from those of refs. [6,19,20] a 
little. However there are some of disagreements between our results and those of 
ref. [21]. 
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Figure Captions 



1. Some Typical Feymann Diagrams for the Subprocess g + g —>■ Bc{B*) + b + c. 
(Here the spots in the figure denote the bound state Bc{B*). In fact, there are 
36 diagrams in total, but here we plot some of them only for illustrating the 
existence of various mechanisms involved in the process.) 

(a) The Feymann Diagrams Correspond to that for Fragmentation mecha- 
nism mainly. 

(b) The Feymann Diagrams Correspond to those for Some else besides the 
Fragmentation Mechanism mainly. 

2. The the Transverse Momentum Pt Distribution for Bc{B*) Productions at 
Tevatron and LHC for Approach-I and Approach-II. 

(a) The Differential Cross Sections of the Bc{B*) Productions versus the 
the Transverse Momentum Pt at Tevatron. The solid line (^4 V) and the 
dashed fine (C V) are those of the B* productions obtained by Approach- 
I and Approach-II respectively; The dashed-dotted line {A P) and the 
dotted line (C P) are those of the Be productions obtained by Approach-I 
and Approach-II respectively. 

(b) The Differential Cross Sections of the Bc{B*) Productions versus the the 
Transverse Momentum Pt at LHC. The notation of the lines is the same 
as Fig. 2 a. 

3. The Differential Cross Sections of the Bc{B*) Productions versus the CM. 
Energies of the Colliding Gluon-Gluon at the Colliders for Approach-I and 
Approach-II. 
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(a) The Differential Cross Sections of the Bc{B*) Production versus the CM. 
Energies of the CoUiding Gluon-Gluon at Tevatron. The notation of the 
hues is the same in Fig. 2a. 

(b) The Differential Cross Sections of the Bc{B*) Production versus the CM. 
Energies of the CoUiding Gluon-Gluon at LHC The notation of the lines 
is the same in Fig. 2a. 

4. The Distributions of the Transverse Momentum Pt and the Rapidity Y for 
Bc{B*) Productions at Precise Energies of the Gluon-Gluon Collisions 

Respectively. 

(a) The Differential Cross Sections of the Bc{B*) Productions versus the 
Transverse Momentum Pt at \/I = SOGeV. The notation of the lines is 
the same in Fig. 2a. 

(b) The Differential Cross Sections of the Bc{B*) Productions versus the 
Rapidity Y at \/I = 30GeV". The notation of the lines is the same in 
Fig.2a. 

(c) The Differential Cross Sections of the Bc{B*) Productions versus the 
Transverse Momentum Pt at = 60Ge\^. The notation of the lines is 
the same in Fig. 2a. 

(d) The Differential Cross Sections of the Bc{B*) Productions versus the 
Rapidity Y at \fi — &QGeV . The notation of the hues is the same in 
Fig.2a. 



18 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408242v3 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408242v3 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408242v3 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408242v3 



